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ABSTRACT 


PetrcOhetrcal investigation of the etfects of short period 
angle of attack oscillation on air-to-ground weapons dive 
delivery accuracy was conducted. The investigation included 
the computer modeling of the longitudinal control system of 
the F-4 aircraft. This longitudinal control system was used 
to fly the aircraft through its standard conventional dive 
bomb delivery maneuver, making corrections for various dive 
angle errors, detected four thousand feet above release 
altitude. 

The dynamics of the F-4 aircraft motion was represented 
by a matrix of dimensional,’ body axis longitudinal stability 
derivatives and force coefficients. This matrix representa- 
tion of the aircraft dynamics was used to solve for the time 
history of aircraft perturbation, velocity, angle of attack, 
pitch angle, pitch rate and flight path angle, following 
COmrection to the desired delivery maneuver flight path. 

The perturbation angle of attack is a measure of the 
gunsigne aimpoint error at any given time, if all other deée- 
livery parameters are met, and converts directly to weapon 
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DEFINITIONS OF SYMBOLS 


REFERENCE QUANTITIES 


FRL 


totaieartreratt velocity in reference state 


aircraft velocity parallel to FRL in 
reference state 


diverareeyelocity perpendicular to FRL 
iierekomence State 


aircraft reference Mach number 


reference flight path angle 
Téetemence pitch angle 
Teference angle of attack 
aircraft mass 

aircraft weight 

air density 

wing area 

wing mean aerodynamics chord 
dynamic pressure 
acceleration of gravity 
moment Of inertia about y axis 
tirust incidence with FRL + oa 


0 


perpendicular distance to thrust line 
inom Cc. 2 


longitudinal reference axis-fuselage 
nome hence lane 


ieiereudinal distance from c.g. to 
bobweight 


(EU Sec) 


(ft/sec) 


(ft/sec) 


(radians) 
(radians) 
(radians) 
(slugs) 
(pounds) 
(slugs/ft*) 
CEE) 

(ft) 

(Cy ata) 
(ft/sec) 
(stuc-f£t.) 


(deg) 


(ft) 


(ft) 





PERTURBATION QUANTITIES 


perturbation flight path angle 
perturbation pitch angle 


perturbation angle of attack 


perturbation velocity component along FRL 


perturbation velocity component perpen- 
dicular to FRL 


perturbation pitch rate 


perturbation stabilator position 


TOTAL QUANTITIES 


ALPHA 


M 


total flight path angle 
ieee y = © - ALPHA 


EQuaieapr1tch angle 9 = 0, + 0 


total FRL angle of attack 
ALPHA = a, + a 


total velocity compert along FRL 
=U, + u 


total Mach number 


(rad or deg) 
(rad) 
(rad/milli- 
radians 
when refer 
to sight 
angle) 


(ft/sec) 


(ft/sec) 
(radians/sec) 


(rad-pos, 
L.ED) 


(rad/deg) 
(rad) 


(rad/mils) 


(ft/sec) 


STABILITY AXIS DIMENSIONLESS COEFFICIENTS AND DERIVATIVES 


D/qS 
9Cy)/ do 


oC, /aM 


9C,/ 95, 





E = L/qs 
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Cy = 8C) /da 

Cio = [2V, /c] aC; / 9 
C = 9C,/0dM 

Lay L 

“La = Le EN OG eR 
Ths = 9C, /96q 

on = m/qSc 

at = 9C_/ da 

oa = [2V_ /c] ac /aa 
CoM = oC / aM 

Cia = [2V_ /c]ac,/3q 
“mS « = 9C_/95o 


BODY AXIS DIMENSIONLESS COEFFICIENTS AND DERIVATIVES 


Cy = Ch cosa, - Cy Sina, 

Cy = (C, - Cy) COSX) - (C), - C, Sind) 
Cyn = CoM COSA - Coy SINQ) 

G = C COSQ,) - C Sind) 

X06 Ddq Lo. 
Cy, = Cy COSQ_ + C)Sindo 

et = (Ce iH Cy) COSA) + (Cy. - Cy) Sina 
G = C COSd) + C SING, 

NO Lé. Dd. 

CAs = C) gSindo 
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CaM = C McOSq =: Cystine, 


All moment coefficients and derivatives are the same for body 


axis as for stability axis. 


BODY AXIS DIMENSIONAL DERIVATIVES 
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CONTROL SYSTEM DIAGRAM TERMS 


pilotestick force 
bobweight stick force 
ci ae 


control stick movement 


stabilator movement due to stick travel 


Sstabilator movement due to stab. aug. sys. 


6 + 6 


SST SSAS 


normal acceleration at bobweight 
perturbation pitch angle commanded 


A — ape 


perturbation dive angle commanded 


iz 


(1/sec) 
(1/sec*) 


(1/sec) 


(pounds) 
(pounds) 
(pounds) 
(inches) 


(radians) 
(radians) 


(radians) 


(£t/sec*) 


(rad) 


(radians) 


(rad) 





I. INTRODUCTION 


The theory and practice of fixed gunsight, air-to-ground 
weaponry is well established and differs only slightly from 
one aircraft to another due to minor handling quality and 
weapons release system differences among aircraft. Just as 
well established is the asa Of “IMpaceeerner Senor, il. 
ties to variation of actual from desired release parameters. 
Table I shows the desired release parameters and error sen- 
Sitivitieswin the most widely used F-4 dive bomb delivery. 
It 18 not the purpose of this study to disagree with any of 
the release parameters for this or any other delivery maneu- 
ver, but rather to demonstrate that increased emphasis in 
Erdininge@ommone Of the error sensitivity categories, that of 
gunsight placement as influenced by FRL angle of attack per- 
Cirmbattonememean gpreatly increase pilot competence in alr-to- 
ground weapons delivery. 

A lack of smshas ig on angle of attack input to gunsight 
Paetimessecems justified on the surface, since everything 
which is an input to sight angle can be computed and set 
prior to the delivery maneuver and should not be subject to 
euroree in addition, Table I shows that the impact error 
flemeowa sight picture error on one mil (0.1% of a radian) 
everett than the impact error due to an airspeed error 
aoe Piots, or a flight path dive angle error-of one de- 


gree. For these reasons and since only one FRL angle of 
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rap rele 


Destred Velivery Parameters and Error 
sensitivities for Maneuver Studied 


Release Parameters 


Dive Angle 40 Degrees 
Release Airspeed * $00 KTAS 
Release Height 5000 Feet AGL 
Sight Angle Om Masks 
Release Horiz Range Py eer ce t 
Release Slant Range 6875 Feet 
Parametet Error: Lupact —Erroc 

Dive Angle 

1° Steep 33 Feet long 

1° Shallow 35 Feet short 
Release Height 

100 Feet High 13 Feet short 

100 Feet Low 1S e@ie Iter 
Release Airspeed 

i@ KTS Fast Soreet Lone 

10 KTS Slow 57 Feet short 
Sight Angle/Aim Point 

+ 1 Mil/l Mil Long fOmreces.long 

- 1 Mil/l Mil Short 10 Feet short 


(This information was calculated from data in Ref. 3. Sight 
angle includes weapon trajectory drop, FRL angle of attack 


and gunsight parallax correction.) 


attack is associated with one steady flight condition at a 
given weight and configuration, it has been accepted that the 
angle of attack input to the sight setting is a quantity which 
varies only with release airspeed and dive angle. This is 
Pe@ermenowever, Only 1f the aircraft is free from all longi- 
tudinal oscillations rising from the most ae ne atmospheric 


PieauGmanece sor 1¢ Longitudinal control action were taken. 
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Most of the emphasis given to sight error due to angle of 
attack perturbations in the past can be summed up in the 
Statement: "The aircraft must be in one-G flight at release 
in order for pipper placement to be an accurate indication 
of bomb impact point."’ Very little emphasis has been placed 
Spe echceudetmmmat a Linite time interval must sexpire atter 
wre ies t "re'ceirt en control “input befere the actual 
FRL angle of attack matches the steady state angle of attack 
for the desired flight condition. If bomb release occurs 
during this finite interval, any difference between perturbed 
FRL angle of attack and angle of attack computed for the 
Sight setting will manifest itself as a bomb impact error. 
Finding the magnitude and duration of this angle of attack 
perturbation after making various dive angle error correc- 


tions was the purpose of this investigation. 
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Pi. PROEEEM SEP -UP 


The standard F-4 visual dive delivery of the MK 82 
LDGP Bomb with conical fins was chosen for the investiga- 
tion. This maneuver involves a roll-in from 12,000 feet at 
400 KTAS on a heading 90° to desired run-in line. Roll-in 
is accomplished by a 135° bank angle descending turn to a 
40° dive along the run-in line. Roll-in is accomplished at 
military power and completed by approximately 10,000 feet. 
Power 15 reduced to idle a@ 430 KTAS (at approximately 9,000 
feet), and the aircraft accelerates at idle power in the 40° 
dive to 500 KTAS at the release altitude of 5,000 feet AGL. 

The largest impact error in range is due to dive angle 
error at release. The causes for this error, shown in 
Figure l, are: 

1. In a dive angle other than 40°, with the gunsight 
pipper on the target, the range to the target is different 
than that planned for at S,000 feet, while the bomb range 
from this point is approximately the same as if in the plan- 
ned dive angle. 

2. The FRL angle of attack is dependent on dive angle, 
but the input to the gunsight for angle of attack is based 
on a 40° dive angle. 

In this study, pilot detection of a dive angle error at 
a0meceot (a typical point to detect an error, shortly 


Pee rerol)-out on final run in line) caused the pilot to 
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make an additional dive angle change to correct for the error 
in dive angle. 

There are two different situations which can be considered 
to be the same dive angle error. A five degree shallow error 
exists for either of the following cases: 

1. The pipper placement is correct for 9,000 feet alti- 
tude but the aircraft dive angle, as noted visually outside 
the aircraft or on the aircraft attitude indicator, is 35°. 

2. The aircraft is actually in a 40° dive but the pipper 
placement at 9,000 feet altitude is 5° short of where it 
should be on the ground at that point. 

In either case the aircraft is approximately 1,850 feet 
(horizontal range) farther from the target than it should be 
at that altitude in the maneuver, and the pilot must shallow 
the dive still further to re-position the aircraft on the 
proper flight path over the ground, and then re-establish 
the desired 40° dive angle. Since pilots rapidly acquire 
Sheed tenolmeimne into a desired Give angic mene normal aiye 
angle error situation is to be in a 40° dive but find the 
pipper aimed too far short or not short enough of the target 
at the 9,000 feet altitude check point. This was the case 


considered in this investigation. 


A. DIFFERENCES BETWEEN ACTUAL MANEUVER AND COMPUTER SIMULATION 
ifeneceulal Maneuver 
Imeran actual run the@aireralt as constantly increase 


airspeed from approximately 430 KTAS at 9,000 feet to 500 KTAS 
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at 5,000 feet with changing lift and drag coefficients. Air 
density increases as altitude ores and eae Career 
erevatiroen 21S not at sea level. 
2. Computer Simulation 

In the computer simulation it was assumed that re- 
lease airspeed was attained at 9,000 feet and no further 
acceleration was achieved during the runs except as a result 
of perturbations due to longitudinal maneuvering. Since the 
assumed quasi-steady state velocity is the actual desired 
release velocity, results in the area of release altitude 
were very accurate. Air density was assumed constant at 
standard day density for 7,000 feet M.S.L., the mid-altitude 
for the run. Target elevation was assumed to be standard 
sea level. Other assumptions and conventions used were: 

a. Flat earth fixed in space, 

b. Atmosphere fixed to earth (no wind), 

c. Rigid airframe, 

d. Body axis system along FRL used throughout, 

e. X-Z plane of symmetry, and 


i> Small perturbation angles and velocities. 


B. REPRESENTATION OF AIRFRAME DYNAMICS 
The following body axis force equation system was ob- 


tained from Ref. 1 in the Laplace demain, 


>) - - X 
S x as cos& XO OMe Md) 2 ESE), u bc 
a sing (UU gtZs)S-2, 7 (Wee vege, a] = “86 
len | 
-M + ~=——T -M> + M Sees Wane M. 
u l. u fo] o} q NI S 


DS, 





where notation and definitions are given in Appendix A. The 
above system is represented in the time domain by the follow- 


ing linear differential equations. 


1. u= OF a T cosé)u + Xa = €COSO ne sage h es “ae +X. 6 





So 5 
: ; U 
; . zy | ie sing Abas ZY _ gsind, te a! : 
Dipoee Zs Uy - 25 U,-45 o, ae 
“6 
+ - 6 
Ue Z* S 
3 Q = q 
° Lin e 
4. qe (MG ti 28 T Ju + M; 17 M = Maa at Noe 5. 


Substitution of Equation 2 for a into equation 4 yields 
the following time domain matrix differential equation in 


the perturbation quantities u, a, @ and q. 








u X tT cosé ee “gEoseo, ae 
. a . =os] U 
‘ Za Tsing Za i Leite. 
mime > 25 Uo 45 Uo 44, Uo 45 
6 0 | 0 0 
: gsind, (Sie 3) 
q M Bache. rea a eit Me 7 -M imewaL MMOD ~~ 
; etl “so 7s ne 0 “a 0 es 
y a 
u aoe 
OS 
a 
URLs 
a 
x + [6 ] 
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aoe 
‘i ee Wee 
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C. COMPUTER REPRESENTATION OF F-4 LONGITUDINAL CONTROL SYSTEM 
Reference 1 contains diagrammatic representation of the 
individual components of the F-4 longitudinal control system 
as a linear control feedback system which can be combined to 
form the complete linear control system shown in Figure 2. 
Either Block Diagram Algebra or the Signal Flow Graph tech- 
nique may be used to reduce the system in Figure 2 to the 
Simplified closed loop system shown in matrix equation form 
in Figure 3. Derivation by Block Diagram Algebra techniques 
of the following Laplace domain representation of the con- 


trol system is presented in detail in Appendix A. 


(1+G2G3G,HiH2(S) - GyGsHs3Hy (S)) 5, (S) = G1G2G3G,(S)[y.-y(S)]. 


As shown in Appendix A, precise representations in the 
Laplace domain of all the control system transfer functions 
would lead to a 42nd order differential equation, the numer- 
ical solution to which would be impractical for the purposes 
of this study. For this reason a method was sought to sim- 
plify the more complicated transfer functions which were not 
directly indicative of airframe dynamics: the aircraft feel 
Syotcmmeaderined Go(S); the stabilator actuator, defined G,(S); 
and the vertical acceleration feedback system to the Bob- 
weight, defined H,(S). These simplifications are shown in 
Appendix A. The 8/5. (S) peamister t£unetion, dehined G- (Sie 
aeemoueh complicated, was not simplified since it represented 
actual airframe dynamic response to stabilator movement and 


its accuracy was considered critical to the investigation of 
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the problem. With the above simplifications, the following 
equation is the Laplace domain representation for the sta- 
bilator position in terms of pilot stick force SG and dive 
angle change required RNG sev Jie (See Appendix A for deriva- 


tion.) 


(S° + 18.458" + 35.57S* + 95.0S* + 1.201S + 0.0224) 6,(S) = 


- 0.01113K, (S° + 3,802S* + 17.99S? + 16.37S* + 1.20158 


0.0224) ly -y(S)]. 


In the time domain this becomes the following fifth 


order linear differential equation: 


O02 Ore 


(S) (4) (3) 5 
6. + 18.456, Booey Oc 2 95.05 — 7 davis F 


5 


= 0.01113K, {(y.-¥) ©? + 3.802(y_-y) + 17.99(y.-y) 


+ 16.37(y.-y) + 1.201(y.*y) + 0.0224 (y,-Y)]. 
iicmaeieemspace technique G£ Reference 5 for representa- 
PIOnmOreinutoOrder, systems of Linear differential equations, 
in which ae forcing function involves derivative terms, was 
used to cast the control system as a set of five linear first 
order differential equations, in general, for the following 


fifth order linear differential equation, 
S 4 3 . ° 
y ds acy! ds acy‘ ds aay = aay * agyY * 
_ ga (5) (4) (3) ; 
B 6 + Bru + Beu + Bou + Bou + Bou 


with the following definition of terms 
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il 
“< 
t 
oOo 
© 
oS 


Xs 

X¢ = X¢ = bou 
Xq = X¢ - beu 
Xe = a = bou 
Xg = ie = bau 


where bg »bc be, br, b, and bg 29@ Gleseijieyal [ny 


Ye By 
be = Be - acby 
be = Be - agbc - agb, 


g ~ agbg - agb, - azbe - agbe - agby 


9 
The matrix representation of the above set of linear 


dafferential equations. is 


5 5 
X6 0 0 1 0 0 X¢ be 
Xo{ = | 0 0 0 1 0 Xo + bo} [ul 
Xo 0 0 0 0 1 Xo be 
Xg | “ag “ag -ay a “ay Xg Po | 


ZS 





XX oX mu + by Reale 


Zot Leo 0 O}{x, 7XgXq 


*6 
iiechemetaemotmrmuencOntrol system, z is defined to be the 


stabilator position 6 and u is the forcing function hee ele 


S? 
The subscripts 5 through 9 were used since the control sys- 
tem equations were to be coupled with the four airframe 


equations of motion to form a ninth-order system for numeri- 


cal solution. 


D. AIRFRAME AND CONTROL SYSTEM COUPLED 

The airframe and control system equation were coupled 
through So; defined by oa on Usine the tact that 
flight path angle [ is defined by If = Q-ALPHA and the flight 
path perturbation angle is defined by y = 8 - a, the finai 
set of nine coupled linear differential equations appears as 


the feedback control system of Figure 3. 
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Pri PROBLEM SOLUTION AND RESULTS 


me CASES CONSTDERED 

Two separate initial cone desione Were trea tedemu nectic 
first case the pilot discovered at 9,000 feet that he was 
2.2 degrees shallow of his desired flight path. In the 
second case he discovered that he 1.9 degrees steep of the 
desired track. For each dive angle error two alternative 
techniques were used to re-establish the desired 40° track 
over the ground. The first was to shallow or steepen 
further an amount which would allow re-establishment of the 
desired track immediately prior to the release altitude of 
5,000 feet. The second technique was to shallow or steepen 
Home thane tmestirst to allow re-establishment on the desired 
MmPcht patneeariy enough for the major protion of the air- 
frame short period oscillations to die out prior to the 
Remeacewartmrude of 5,000 feet. Figures 4 through 7 illus- 
Miitemellcmiwocrror cases and the flight paths used to cor- 


rect for these errors. 


B. METHOD OF AIRCRAFT CONTROL 

In shallowing or steepening the aircraft dive angle, 
the full system of nine equations was used to solve for the 
time response of u, a, 8, q and y to the commanded perturba- 
jeon dive angle liar When the correction dive angle or the 


re-established 40° dive angle was approached, the reduced 


tai 
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System consisting of only the four airframe equations was 
used to solve for the perturbation quantities. This was 
done for two reasons: 

1. The simplified transfer functions of the longitudi- 
nal control system were invalid in the pitch rate range be- 
tween -0.5° per second and 0.5° per second. Obviously a 
constant dive angle or -zero -piteh-rate falls-wrthin this 
range. 

Z. Simee it was the short period angle of attack per- 
turbation which was of interest, the radian frequency for 
this mode of oscillation was such that airframe damping 
mBavhner thal pilot control movement 1S the factor which 
dampens the oscillation. - 

Theguatoetecnen flew the aircraft from the initial dive 
angle to the correction dive angle through the control sys- 
tem, anticipated achieving the new desired dive angle based 
on pitch rate and set the stabilator to the reference posi- 
tion (perturbation S¢ = 0) for the new dive angle. 

After the aircraft had flown the required time on the 
adjusted flight path to approach the correct 40° dive angle, 
the controlling sequence was repeated with the pilot re- 
establishing the 40° dive and the airframe dynamics in the 


new reference state maintaining it. 


C. COMPUTER THEORY FOR AIRCRAFT CONTROL : 
Reference 6 contains computer programs for computational 


assistance in the study of linear control theory. The 
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programs used were GTRESP, CALCU, RUNGE, TRESP, YDOT, Y8VSX 
and user supplied subroutine RFIND. Program listings and 
explanation are contained in Appendix D. The purpose of the 
programs 1s to furnish the graphical time response to the 


following system: 


{x(t)} = [A] {x(t)} + {b} [u(t)] 
fu(t)] = K(r(t) - tk}! ec Gao 
y(t) = {c}? {x(t)}. 


In the case under consideration, {x} is the vector 
of perturbation quantities, [A] is the nine-by-nine or 
four-by-four plant matrix defining the airframe and control 
system, {b} is the control vector, K = K is the pilot stick 
force per radian of difference between actual and commanded 
dive angle, r is the commanded perturbation dive angle, {k} 
is the vector of feedback coefficients, and y is the actual 
perturbation dive angle made up of the output vector {c} 


multiplied by the airframe perturbation vector 
zs T 
(y = {0 -1 1 O0}{u a 6 q} ). 


Following is the detailed breakdown of one of the 
Situations considered. At time zero the pilot at 9,000 
feet altitude in a 40° dive determined that he was 885 feet 
horizontal range short of the correct position for that al- 
titude in the delivery maneuver. This was equivalent to 
being 2.2° shallow. He elected to correct for this by shal- | 
iewane 5°, flying at a 35° dive angle until approaching the 
desired track, steepening back to a 40° dive angle and flying 


meme 20 angle to release. 
Os 





GTRESP was called with the following data: 

1. the nine-by-nine plant matrix with airframe data 
referenced to a 40° dive angle and the control system, 

7 eeme -Conenrol vector, 

3. the output vector, 

4. the feedback coefficient vector, 

Se SG based on required dive angle change of 5°, 

6. ve equal to +5° in radian measure, and 

7. ezeroeimitial perturbation vector referenced to 
40° dive conditions. 

GTRESP then used the other sub routines tO pPperrorned 
fourth order Runge-Kutta integration to compute and plot the 
firme YTesponsesore the Nine perturbation quantities. 

At 0.975 seconds the perturbation dive angle was 
approaching 5° or I was approaching 35°. The airframe per- 
turbation quantities for time equal to 0.975 seconds were 
transformed to perturbation quantities for a 35° dive ref- 
erence state, and GTRESP was called with the following in- 
put data: 

1. the reduced plant matrix of airframe data only 
referenced to a 35° dive, 

Zee OUtpuUt VYECtor, and 

Mecirerdtteperturbatrion qUuancities referenced £0 
35° dive at time equal to 0.975 as the initial condition 
Viet On. 

Poieinne to input the control vector, feedback coeff1- 


Cake mies. - and ie bypasses the control system which 1s not 
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accurate for zero and very small pitch rates and simulates 
the pilot controlling the aircraft ee and maintain- 
ing the zero perturbation stabilator position for the 35° 
dive reference condition. Again GTRESP computes and plots 
the time response of the aircraft perturbation quantities. 

At 7.2 seconds total elapsed time the aircraft was 
approaching the desired’40° track over the ground and GTRESP 
was called with the large plant matrix for control back to 
a 40° dive. Input was the same as for the first maneuver 
except that the aircraft dynamics portion of the nine-by- 
nine matrix was based on the 35° dive since that was the 
condition from which perturbations occurred, and the ini- 
tial condition perturbations remained referenced to the 35° 
Gaiay.e 

GTRESP computed and plotted the perturbation quan- 
tities through the transition back to a 40° dive. At 8.175 
seconds total elapsed time the 40° dive on the correct 
flight path was re-established and it was simulated that the 
pilot neutralized longitudinal control by calling GTRESP, 
with the reduced plant matrix and initial perturbations 
transformed to those for a 40° dive reference state, to 
ealculatesand plot the perturbation dive angle, velocity, 
mieeomoneattack, Pitch angle, and pitch rate to release 
Ode, . 

GOnethis particular run it was arranged to have the 
@ineratt re-achieve the desired 40° dive angle at exactly 


Pre 5,000 feet release altitude. Therefore, the angle of 
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attack perturbation, converted to mils, existing at 8.175 
seconds is in fact the release sight picture error that the 
pilot would have had on this maneuver at release altitude. 
Pertinent aircraft position and gunsight error information 
concerning each of the four bomb delivery maneuvers simu- 
lated is listed in Table II. 

Computer output showing time history in tabular form 
of aircraft perturbation during the four phases of each of 
MieecOUr Tunsetested and in graphical form for the final 


delivery phase is enclosed in the Computer Output Section. 


D. RESULTS OF THE FOUR RUNS 

Table II shows that the total time available for ma- 
femavering @etrer detecting a dive angle ecumoen, 4,000 feet 
above release altitude, is very small for a 40° 500-knot 
dive bombing run. Therefore, any corrections made must be 
both timely and correct, if they are to succeed in increas- 
ing the accuracy of the bomb impact. Table II also shows 
ideo vemenreiy increasing the correction dive angle from 
5° to 8°, and thereby re-establishing the desired dive path 
1.75 seconds prior to release rather than at release for 
the shallow run and 1.4 seconds prior to release rather than 
at release in the steep run, reduces the gunsight position 
error from 83 to 9 mils for the shallow case and from 80 to 
omits tor the steep case. 

Peamination of the graphical time histories of perturba- 


tion angle of attack, presented in Appendix C, shows that 


36 








"2IV OSseoToy 3e LO1Ig YYstTSsuND pue SUuOT}RI9IIOD) oTSuy oATG BuTIing uotiIsog WFeIDIITY “TI STqey 

¢ §26°9 0 S78 0 Coe OPTTO .Or 

- GZS°S CL8 002 Coe L¢o9 Para Oo) alice Sash 

= SLE°P BIE SM C8Zt 69¢CT O86T Sp Pl oeeou 

= SZc f 9g¢¢ SZL 67c¢ TS9 o8t OF TSAO YSng 

uOT}399I11049 as 

d903S ,6°T 

b NNY 

08 = - = : z F9PITO Or 

- SZ°9 0 019 0 6SS o0v 02 dn [Ing 

- SLL 5 0OT9 C982 6SS PS88Z oP Fear 

= S76. CLVE 809 Crrse LSS oS OF TSAO YSN 

UOT] Dor 0) 2G 

do02S .6°T 

Cc NANG 

6 S22 °8 0 C2GE 0 0S6 PPTL) OP 

= SZS°9 CELL 6CL 0S6 0729 o0P OF ADAO YSN 

2 c¢o’°s TZ8T 9682 OLST 681 SPliomaec. 

= Soc EL LOLD v88 68¢°¢ TT9 E7S Os Gne pia 

a) sales) 3c Ica 08 

MOTTBUS ,2°2 

¢ NN 

cs = - = = = 2PTTS 30 

< SZT°8 0 $S9 0 voP o0V OF TOAD YSN 

= COL SS9 VVCP v6P~ 900¢ IPL LO eas 

= SZ6° 666? cS9 O00S¢ 00S oof 0% dn [Ing 

(sT Tw) (99S) (at) (23) Ce) es) UOT}901109 85 

OSBOTOY SUT] OSBOTOY 02 POTToOAel], SSBOTSOY S07] MOTTEUS sae 

1B posderTy 9due.stg adueYSsSTG dAOQGY an vf T NOY 
Onlee] TRIO] TeUOZTIOY TeUOZTIOY VYUSTOY 


on 





any additional time spent on the desired dive path rapidly 
decreases the angle of attack perturbation, thereby rapidly 
increasing the accuracy of the gunsight position relative 


to actual bomb impact point. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


Table II shows the elapsed time in each of the four 
weapons delivery maneuvers tested at which the aircraft 
finally becomes established in a 40° dive on the proper 
track over the ground for an accurate bomb hit. The Com- 
puter Output shows the FRL perturbation angle of attack at 
numerous time increments in each run. It is obvious for 
both the initially steep and initially shallow runs that 
the 5° correction dive angle, although sufficient to correct 
mie, adircrart back to desired dive path by release altitude, 
1s not large enough to correct back to the desired dive 
path soon enough to allow for sight picture settling time, 
due to angle of attack oscillations of the short period 
longitudinal mode. However, by merely increasing slightly 
the magnitude of the correction dive angle change, it was 
possible to reposition the aircraft on the desired flight 
path early enough before the release point to allow time for 
Miewsight picture to stabilize sufficiently to achieve an 
acceptable bomb impact. Examination of the graphs showing 
time history of angle of attack perturbations during the 
final phase of the delivery maneuver shows that even earlier 
establishment on desired dive path is desirable and would 
increase the accuracy of the sight picture further. 

fmectetual practice, a pilot on detecting a dive angle 


error would not consciously attempt to hold a dive angle 


he 





correction for a specific amount of time with the idea that 
imeso doimlgeie would re-establish his aircraft on the proper 
flight path. If he could do this, there would be no need 
for the gunsight. All the cases described herein would 
actually achieve very accurate bomb impacts, since for an 
aircraft at the correct velocity and in the desired dive 
angle along the proper path over the ground, very large 
angle of attack perturbations produce bomb impact errors 
only insofar as they change the aircraft dive angle. The 
error which would actually be produced in the case where 
all delivery parameters are met except for an angle of at- 
tack error would be that the gunsight reticle would be far 
Ont the tanger. 

However, since the pilot in a visual dive bomb delivery 
Mees Visctialiyeacquired data for dive angle and aircraft 
poonetonineesne 1S Entirely dependent on his gunsight for 
accurate aircraft positioning information; and during a 
period when he is receiving inaccurate gunsight information, 
he has no good indication of his position relative to the 
Geeared release position. 

It has long been standard practice in visual air-to-air 
gunnery with a gyro-stabilized gunsight, which develops its 
own lead angle based on the normal acceleration of the fir- 
Mem isenatt required to track the target, to teach students 
fmoemcomsaueeze the trigger until they have maintained a 


emoormetracking solution for two seconds on the target. This 
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= to allow angle of attack oscillations from a changing 
pitch rate to die out prior to firing. 

The need to apply this same technique to visual air-to- 
ground delivery is just as important, or even more so when 
it is considered that a gyro-stabilized sight has inertia 
to decrease the magnitude of the effect' of aircraft angle 
of attack perturbations; whereas the sight line of a mech- 
anically fixed sight moves directly with angle of attack 
perturbations an amount equal to the movement of the FRL. 
Nor does a pilot in air-to-ground delivery have the option 
to wait until perturbations die out if this will take him 
below the planned release altitude, since the computed sight 
angle 1s good for only the planned delivery parameters. 

What he must do 1s to detect dive angle errors sufficiently 
early in the run and make corrections sufficiently large 
that all corrections are completed and all significant 
Pertrlnvattonsscliminated prior to the planned@release hereht. 
iticmeriiay spentained specifreally togthe F-4 aircraft in 
a 40° dive, MK 82 LDGP Bomb, delivery. However, since all 
Navy eiwehter and light attack jet aircraft have similar 
short period characteristics with slight variance in damping 
Rartowamae oscillation frequency, and since all fixed-sight 
air-to-ground weapons delivery maneuvers have similar error 
Seucrervities (differing only in the magnitude of the range 
Paaometor cach mil of aimpoint error), the results are ap- 
plicable to all fixed gunsight, visual, air-to-ground weapons 


@eliveries in all Navy tactical aircraft. As such, it 1s 
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recommended that similar studies be made to encompass all 
air-to-ground weapons deliveries from Navy tactical aircraft 
and the results be used in the lecture portions of air-to- 
ground training in both fleet and replacement squadrons. ihe 
addition a graphical display of the results herein can be 
easily demonstrated to a pilot in flight in the following 
Manner. Trim the aircraft for Steaty"revel “flicht With the 
gunsight set on the distant horizon. Make a small but sharp 
longitudinal displacement of the control stick and return it 
to neutral. Notice that, although the change in aircraft 
track (speed and altitude) is extremely small, the magnitude 
of the gunsight oscillation above and below the horizon is 
Significant; and that while the time to damp is short, it is 
Hong enough to lose several hundred feet, if it were in a 
40° S500-knot dive. 

Awareness of this phenomenon can do much to improve the 
proficiency of the not atypical pilot who, on being ques- 
tioned about pipper placement after an unsuccessful air-to- 
ground weapons training flight, is not really able to tell 
an instructor whe-<c his pipper was positioned at release on 
most of his runs, this because of late corrections resulting 


in receipt of meaningless data from the gunsight. 
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APPENDIX A: DERIVATION OF LINEAR CONTROL FEEDBACK COMPUTER 
REPRESENTATION OF F-4 LONGITUDINAL CONTROL 
SYSTEM 


The longitudinal control system of the F-4 aircraft is 


Shown below in diagrammatical form, 


= 


on ee = pea ee ee ee ae 





where, 
oo) Syecen = G, (S)eay a = 5) ee 
Sato S560 5.5 poun orc 


= __ radians stabilator travel 
Control Stick Gearing = G,(S)= 0.05864 _ ({ mere ERREET CE ¢ TAVel DM 


It 


Stabilator Actuator = G,(S)= 20 [ - 3 ] 


Bobweight 2 H, (S)= 0.1663 ( pounds ag Ou eI En force ) 


2 Rad2ans. stabil ateonetwaved 


_ S 
Seab. Aug. System = H,(S)= 0.15 [ =>] ( yadian/sec. pitch rate 


With the addition of airframe dynamics and a pilot, modeled 
by Ky pounds of stick force per radian of error between actual 


ma@emdesired pitch angle, the representation becomes 
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where the additional transfer functions are, 


Pilot force/radian eg = G,(S) = K (pounds force) 
P radian 
ae/6 transfer function 
oe ss 78S +22 72540 04550" — fe 
7 H, (8) = 861.5 4 3 2 maa eae 
St oom bolo ome ee 7 OSt) 50224 


fo, transter function 


S 
2 ; 
= G.(S) = -24.64 [—p———— eT racian) 
St 2 oui omer 16 tO Sr tl. 17 9S+0.0224 
; _ 5 radian/sec 
q/®@ transfer function = H,(S) a ( FSET ) 


- A block diagram algebra reduction of the control system 


Pemococmeped in Reference 5 follows: 








ae er ewww ee > 
A a 









= Je 


HH 


£46,636, 





6 


eee? CS 


= G,G,G,E, 














— eee 0 ee oe 
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G,G 
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8465 
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Sg 2G 


1+G,G6,G,HjH, 





., SsHsH, 
st * “G6; 8s 
= G5G63G, 
1 jdt G,G,G,H,H, 
G5G5G, 
T+ G56,6,H,H, 
a 650,0,6,1.A 
: > G50, (1*G,6,6,H,H,) 
G5G2G, 
I+ G5G,6,8,1, - 0,0,HF, 





GG iid 


Pans 3.44 





me ree pl 














Geen 
2 45 Ded. 


S oT 1 + G,G,G,HjH, = G,GcH ZH, 


G,6,6,6,(@ - 8) 


S 1s G,G,G,H,H,- G,G HH, 


cic Cr Ganas- G,G.HH,) o.(S) = Gigae a, | (cman 


In this final form the control system in the Laplace 
domain can be readily transformed to its state space repre- 
sentation of a linear differential equation for: 5. as a 
function of time and (6. = 9). Howe Ve pe redteiing salieetnc 
transfer functions which are polynomial fractions in S to a 
least common denominator would result in a forty second 
order linear darterential Equation. Therefore, representa- 
tion of some of the transfer functions that would yield a 
resulting differential equation more amenable to numerical 


solution but still accurate, was sought. 


7 (Sj G_(S) was not disturbed since it represented 
Bae airframe dynamics which was under investigation. 
Neither was the Stab. Aug. System = H,(S). Howe Vie bd eee = 
Hieweyeresponse investigation of the control stick = G,(S), 
Pmeestabiilator actuator = G,(S) and ae 


Z 
Peeor Computer program of Reference 6, yielded Bode Diagrams 


= H,(S), using the 


of transfer function magnitudes and phase angles showing 


eiateover the pitch rate range of interest (0.01 to 0.2} 
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radians/second - a normal acceleration range up to + 5.2) the 
above transfer functions could accurately be represented by 
constant magnitudes with zero phase angle. 

With these simplifications, the transfer functions are 


as follows: 


6, (8) = K, H,(S) = 1600 
G,(S) = 0.04792 H,(S) = 0.1663 
G,(S) = 0.05864 H,(S) = S 
_ 7 1 
SCS) =o H,(S) = 0.15[ = 77 | 
G2 1 a5 yay, 
G.(S) = -24.64 [ - - 


Wim a ase... ee” rn era ane: 
Wil es0oGe 4 15, 19S + 1.1799 2 0.0274 


The Laplace domain representation for So as a function 


of time and ce - 6) becomes 


4 5 2 


fovesamicmensic’ = 35.575" + 95.0S 


4 


Zee Ots =) 00224) 6. (S) = 


5 Z 


: 0.01113K, (S° + 3.802S° + 17.99S”° + 16.37S° + 1.2015 


mmoeoe24)((0. - 6(S)]. 


item ite Ot Ka? Pulot Loree per radian of pitch angle 
.error, was varied depending on what gain was required to 
achieve an average pitch rate of 0.115 radians per second 
which is equivalent to a change in normal acceleration from 
steady state conditions of 3g (typical in fact of the pitch 


rate a pilot would use in making a dive angle correction). 
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In addition, the pilot was assumed to be aware that 
actual flight path angle (dive angle) lagged the pitch at- 
titude angle during a changing pitch rate by the amount of 
instantaneous angle of attack. Since the purpose was dive 
angle change not just pitch attitude change, 6-a or y was 
fed back in the control system to solve for 6. as a function 
of time and gah rather than (6-8). This is the repre- 
sentation shown in the matrix feedback differential equation 
of Figure 3. 

The simplification of the transfer functions listed 
above necessitated using airframe dynamics alone, in a 
Simulated stick fixed condition, rather than coupled with 
the control system after the desired dive angle was approach- 
ed, since very small pitch rates were outside the limits for 


which the simplified forms of G4(S), G,(S) and H,(S) applied. 
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APPENDIX B:; NUMERICAL VALUES OF CONTROL VECTOR AND ELEMENTS 
OF THE PLANT MATRIX FOR REFERENCE DIVE ANGLES 
USED 

The elléments of the control vector derived in Sections 
II-B and II-C are invariant with dive angle. Values of 
elements of the control vector are listed in Figure 8. 

Derivation of formulae for the elements of the plant 
Matrix is in Sections II-B and II-C, and the matrix is 
Shown in symbol form in Figure 3. Definitions of the ele- 
ments is in "Definitions of Symbols" at the front of this 
study. Calculation of the airframe dynamics entries was 
accomplished in the following manner: 

1. Dimensionless stability axis data were obtained 
trom References 1 and 2Z. 

2. Dimensionless body axis figures were computed from 
the "Definitions of Symbols", using stability axis data 
and reference state angle of attack. 

s Dimensional body axis figures were calculated using 
dimensionless body axis figures, dimensional derivative 
iiinrerons, and applicable reference state data. Tables 
III and IV are compilations of the unperturbed 40° reference 
state conditions, and aerodynamic coefficients/stability 
derivatives for each reference dive angle used. Figure 8 
lists the values of the plant matrix elements for each dive 


angle. 
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Since the aircraft was assumed in a constant speed 
(except for perturbations) 500 knot dive, rather than ac- 
celerating from 430 KTAS to 500 KTAS, lift and drag coef- 
ficients from References 1 and 2 were invalid. These 


coefficients were calculated for each dive angle from free 


body force balance diagrams. 
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Table III. 40° Dive Reference Conditions and Stability Axis 
Aerodynamics Coefficients and Stability Derivatives 


To . = -40° = -0.6981 radians 
Ver = 844.4 ft/sec 
0 
U, = 844.4 ft/sec 
ie = 52509 £t/sec 
M, = .775 
or = -39.77° = -0.6941 radians 
Ot = 3.99 mils = 0.2286° = 0.00399 radians 
m = 1210 slugs 
W =. 38920 1b 
0 = 0.001933 slugs/ft” 
S = §30 ft’ 
e = 16.04 ft 
aq 9S ae). i Sey eae 
g apa San Be ft/secé 
I, = 122200 slug fr? 
ak Snes 
Lin = -0.57 ft 
1, = 39.4 ft 
Ch = 0.06535 
Che = 0.1464 
Cou = 0 
“Dé = 0 
Cy = 0.1391 


Sl 
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tives. 
r -32° -35° 
Cy 0.05026 0.0611 
Cy -0.002913 -0.002064 
a 
Cyn 0 0 
6 -0.001416 -0.001416 
Xb, 
Cy 0.1259 0.1317 
Ne 3.411 Sac 
Cys Oey (Ee 
Com 0 0 
Na DZ etee PROD 
ONS. 0.3555 05555 
x, -0.03541 -0.04345 
XO Oy Sel) 0.4344 
xo -~0.0254 -0.0254 
XO, 7 4 On4274 
Zh -0.08029 -0.08759 
Z eee =» 1022. 
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Z° ~2.007 -2.007 
OF 
He -~6.365 -6.365 
25 -107.2 -107.2 
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M -13.76 -13.76 
OF 
Mes -0.4553 -0.4553 
Ne -1.079 -1,.079 
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Body Axis Force Coefficients and Stability Deriva- 


-45° SANE 
- 0.07144 O77 7 
-0.006422 -0.01468 
0 0 
-0.001416 -0.001416 
0.1495 0.1628 
3.432 3.437 
0 0.7 
0 0 
DE. Dee 
0.355 0.355 
-0.05129 -0.05566 
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~0.0254 ~0.0254 
0.4274 0.4274 
-0.1061 -0.118 
- 1045 - 1056. 
52 (Ny 22007 
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109) 2 =i). 2 
SOMO Ores 7 ORO 0d sale 
-13.76 ee 76 
-0.4553 a) dl 35 
=a OS -1.079 
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APPLICABLE TG ALL DIVE ANGLES 


ELEMENTS OF THE CONTROL VECTOR, 
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8 DEGREE CORRECTION 


2-2 DEGREES SHALLOW, 


RUN 2, 


TIME RESPONSE QF PERTURBATION QUANTITIES DURING PULL UP FROM 40 TQ 32 DEGREES 
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TIME RESPONSE OF PERTURBATION QUANTITIES CURING 32 DEGREE OIVE 
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TIME RESPONSE OF PERTURBATION QUANTITIES DURING PUSH OVER FROM 32 TO 40 DEGREES 
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TIME RESPONSE OF PERTURBATION QUANTITIES DURING 40 DEGREE DIVE 
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TIME RESPONSE OF PERTURBATION QUANTITIES CURING 40 DEGREE DIVE 
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TIME RESPONSE OF PERTURBATION QUANTITIES CURING 48 DEGREE DIVE 
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TIME RESPONSE OF PERTURBATION QUANTITIES CURING 40 DEGREE DIVE 
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